EXPECTED LENGTH OF A PRODUCT OF RANDOM 

REFLECTIONS 



JONAS SJOSTRAND 

Abstract. We present a simple formula for the expected number of inversions in 
a permutation of size n obtained by applying t random (not necessarily adjacent) 
transpositions to the identity permutation. More general, for any finite irreducible 
Coxeter group belonging to one of the infinite families (type A, B, D, and I), an 
exact expression is obtained for the expected length of a product of t random 
reflections. 



1. Introduction 

In 2000, Eriksson et. al. [6] studied the expected number of inversions after t random 
adjacent transpositions applied to the identity permutation. Their main motivation 
came from the area of molecular evolution where a genome (modelled as a permu- 
tation) undergoes random mutations (modelled as adjacent transpositions). Since 
then several people have analysed this and similar problems, motivated both by 
biology and pure mathematics. 

In 2005, Eriksen [8] improved the results from [6] and obtained an exact formula 
for the expected number of inversions. More recently, in 2010, Bousquet-Melou [2] 
gave a completely different exact expression for the same thing! In 2002, Troili [7] 
generalized the question to other groups than the symmetric group, and she gave 
corresponding exact formulas for the dihedral groups ^(ui). 

In 2004, Eriksen and Hultman [1] modelled the genome rearrangement problem 
closer to reality by considering any random transpositions (and not just adjacent 
ones) as mutations. Instead of the number of inversions they computed the expected 
absolute length of the resulting permutation, i.e. the minimal number of transposi- 
tions whose product is the permutation. In a sequel [5] they generalized their results 
to a class of complex reflection groups. 

For an expose of the research on more general aspects of Markov chains generated 
by random (adjacent) transpositions, such as mixing times, we refer to Bousquet- 
Melou 0. 

In the present paper, we will generate permutations in the same way as Eriksen 
and Hultman, i.e. by random transpositions, but we will meausure them by the 
number of inversions as in the original papers. We will also generalize the framework 
to all finite Coxeter groups and present exact expectancy formulas for the groups 
A n , B n , D n and I^irn). In Coxeter group lingo, we compute the expected length of 
a product of t random reflections. 

The paper is organized as follows. After a section with some basic definitions, 
we review the previously known results mentioned above so that the reader gets a 
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feeling for what kind of expressions may appear for the expected (absolute) length 
and also what type of combinatorics is going on behind the formulas. Next, in 
section H] we loosen up by computing the expected length and absolute length of 
a product of random reflections or simple reflections in the dihedral groups ^(m). 
After this relatively simple task, in section we are ready to state and prove our 
main results: formulas for the expected length of a product of random reflections in 
the Weyl groups of type A, B, and D. Finally, in section[U]we summarize the current 
state of affairs and reveal some of our ideas for future research. 



2. Definitions and notation 

We will follow the notation from the book by Bjorner and Brenti [1] where also the 
definitions of the basic concepts may be found. For the reader's convenience we have 
included the most relevant definitions below. 

A Coxeter group is a group W with a presentation of the form 

(S | (ss') m(s ' s,) = id for all s,s' e S with m(s, s') ^ oo), 

where m is a function from S to {1, 2, . . . , oo} such that m(s, s') = m(s' , s) for all 
s,s' E S and m(s, s') = 1 if and only if s = s'. The elements of S (which we identify 
with their image in W through the presentation) are called generators or simple 
reflections and the pair (W, S) is called a Coxeter system. The group elements that 
are conjugate to some simple reflection are called reflections and the set of reflections 
is usually denoted by T. 

By definition, any group element w € W is a (finite) product of simple reflections, 
and the minimal £ such that w is a product of £ simple reflections is called the length 
of w and is denoted by £[w). Similarly, the minimal £' such that w is a product of 
£' reflections is called the absolute length of w and is denoted by £'(w). 

Now, let us make a very general definition. 

Definition 2.1. Given a (multiplicative) monoid M , a nonempty finite subset R C 
M , a "length" function <p from M to N = {0, 1, . . . }, and an integer t > 0, we define 

r 1 ,...,r t eR 

to be the expected "length" of a product oft random elements, chosen independently 
and uniformly from R. 

As a special case, for a Coxeter system (W, S) with S finite, E)Y e (t) denotes the 
expected length of a product of t simple reflections. Similarly, if W is finite, E^ e {f) 
denotes the expected absolute length of a product of t reflections. 

As mentioned in the introduction, exact formulas for Egg(t) have been found 
for the symmetric groups (W, S) = A n by Eriksen [3] and for the dihedral groups 
(W, S) = l2{m) by Troili [7]. Exact formulas for Ejfg,(i) have been found by Eriksen 
and Hultman, first for the symmetric groups [1] and later for the complex reflection 
groups G(r, l,n) where the hyperoctahedral groups B n are included. 

In the present paper we will compute for the three infinite families of finite 
irreducible Coxeter systems A n , B n , D n , and I^im). 
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3. Known results 

To get a flavour of the kind of expressions that may appear, we will begin by- 
presenting the already known results and discussing briefly what methods were used 
to obtain them, before we reveal our own results in the next section. 

Theorem 3.1 (Eriksen (2005), Bousquet-Melou (2010)). Let n > 1 and t > 0. The 

expected number of inversions after t random adjacent transpositions applied to the 
identity permutation in A n is 

^<<>-e;K:)e(::;)<-i>~<~*-. 

r=l x ' s=l x ' 

where 

An alternative formula is 

* . . n(n+l) 1 (cos a, • + cos at) 2 / 4 . .V 

k,y—0 

where a k = (2k + l)ir/(2n + 2). 

The proof of the first formula, due to Eriksen, is built on a discrete heat flow 
process in a two-dimensional lattice. His expression has a combinatorial interpre- 
tation in terms of lattice walks. The second formula, due to Bousquet-Melou, was 
obtained by attacking the same heat flow process by generating functions and the 
kernel method. It has the advantage that the dependence on t is very explicit which 
makes it possible to analyse what happens in the limit when n and t tend to infinity 
at the same time in various ways. 

Theorem 3.2 (Troili, 2002). Let m G {2,3,... ,00} and let t > 0. Then the expected 
length of a product oft random simple reflections in the dihedral group ^{fn) is 



where 




Here, binomial coefficients with b negative (or b = —00) vanish by definition. 

Troili's proof is built on the observation that 

t-i 

E%} m \t) = ^[Prob^(TrM) = 0) - Prob(f (vr^) = m)], 

r=0 

where ir^ is a product of r random reflections in Izirn). 
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Theorem 3.3 (Eriksen, Hultman, 2005). Let r,n be positive integers not both equal 
to one. Then the expected absolute length of a product oft random reflections in the 
complex reflection group G(r, l,n) is 



1 ^ lmin ^-^ / r ((P) + _ («-V +2 ) + n) 




n 



>(©+®-rv fl )+p 



where 



and 



(-1) 



6, 



n-p- q +i (P ~ g + l) 2 (n\(n-p-\ 



(n — q + l) 2 (n — p) \p/ \ g — 1 
(_l)n-p-g+i / n \ /n — p — 1 N 



n — p \p J \ q— i 
Note that A n _i ~ G(l,l,n) and £ n ~ G(2,l,n). 

The proof of Eriksen and Hultman's formula involves group representations and 
characters. To see why this is natural, think of the symmetric group ^4 n _i where 
the absolute length of a permutation is n minus the number of cycles in the cycle 
representation of the permutation. Applying a random transposition either splits a 
cycle or merges two cycles, and it is easy to see that during this process we only 
have to keep track of the cycle type of the permutation, i.e. its conjugacy class. 

4. The dihedral groups 

As we have seen, computing E 1 ^^ {t) for the dihedral groups lead to a rather 
cumbersome expression in Theorem 13.21 In this section, we will see that much nicer 
formulas exist for E^\t), E^\t), and E%fr\t). 

Theorem 4.1. For m > 2 and t > 1, the expected length E^ m \t) of a product of 
t random reflections in the dihedral group I2 (m) is m /2 (for all positive t) if m is 
even and y — 2m */ m ^ s odd- 
Proof. Note that in ^(m) the set T of reflections is precisely the set of elements 
with odd length. Thus, multiplying any element in T by a uniformly chosen random 
reflection yields a uniformly chosen random element in l2{m) \ T, and, similarly, 
multiplying any element in I2 (m) \ T by a uniformly chosen random reflection yields 
a uniformly chosen random element in T. A simple calculation reveals that the 
average length of an element in T is m/2 if m is even and y + ^- if m is odd, and 
that the average length of an element in l2(m) \ T is m/2 if m is even and y — 2m 
if m is odd. □ 



Theorem 4.2. For m G {2, 3, . . . , 00} and t > 0, the expected absolute length 
E%p(t) of a product of t random simple reflections in the dihedral group hirn) is 
1 if t is odd and 

2'- 1 ^ \\-km 



\k\<\tl2m\ 
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if t is even. (If m = oo, the sum contains only one term.) 

Proof. Clearly, in a dihedral group, the absolute length of an element w is 1 if £(w) 
is odd and 2 if i(w) is even, unless w = id. The probability that a product of an 
even number t of simple reflections is the identity element is precisely the sum above 
divided by 2*. □ 

Theorem 4.3. For m > 2 and t > 1, the expected absolute length E^y(t) of a 

product oft random reflections in the dihedral group ^('m) is 1 if t is odd and 2 — — 
ift is even. 

Proof. The theorem follows from the same reasoning as in the proof of Theorem 14.11 
and the observation that the absolute length of any element in T is 1 and the average 
absolute length of an element in I2 (m) \ T is 2 — ^- . □ 



5. The Weyl groups of type A, B and D 

In this section we will obtain exact expressions for EjfJt) when W is a Weyl group 
of type A, B, and D. The formulas and their proofs turn out to be quite similar for 
these three types, but the symmetric groups happen to be a bit less technical to 
deal with, so we start with those. 

5.1. The symmetric groups. 

Theorem 5.1. For n > 2 and t > 0, let ir® E A n —\ be the random permutation 
obtained by applying t random transpositions to the identity permutation. Then the 
following holds. 

(a) The expected number of inversions in is 

n(n-l) (n + l)(n-l)/ 2 V (n-l)(n-2) 



T ' e w 4 6 \ n-lj 12 V n - 1 

(b) For 1 < i < j < n, the probability that irf^ > tTj is an inversion is 
1 7 — i ( \ 2 \ * f j — i 1 



2 n \ n — 1 / V n 2/V n 



Proof. Part (a) of the theorem follows from part (b) by the simple observation that 
Ei<> t<j <n(j-i) = n(n 2 -l)/S. 

It remains to prove part (b). For every 1 < i,j < n and t > 0, let 



pg := Prob(*?> < *f) and ug := Q $ 



t 



Clearly, ufj = and for % ^ j, 



Now, define vfj := ug —uf] for 1 < i,j < n. The above recurrence for u translates 
to the following recurrence for v. 



■C- ;>S-^+E(^-^)+E(»S--»S: 
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Using that v, 
(2) 



(t) 



for any i, we can write the recurrence in a more elegant way: 



(t+i) _ n ( n ~ 5) (t) 



i'=i j'=i 

Next, we will express the above recurrence in terms of linear operators. Let A n 
denote the vector space over M of all real antisymmetric n x n matrices, and define 
the linear operator Q on A n as follows. For any v € A n and for any 1 < i,j < n, 
let each entry in Qv be its row sum plus its and column sum in v, or formally, 

n n 

i'=l j'=l 
Now, the recurrence relation ([2]) can be written 

= (Q + X l)v®, 



where I is the identity operator on A n and x = ra ^" 2 5 ~ ) . The starting point v t - [)> is 
the n x n matrix with zeroes on the main diagonal, ones above it and minus ones 
below it. From the simple but crucial observation that Q 2 = nQ we deduce that 



,(o) 



,(*) 



t„(o) 



Since {Qv 



(oh 



ri,j 



(Q + xlfv 
2{j — i) we obtain 



{ n + x yQ +x t (j_Q 

n \ n 



,(0) 



(*) 
v- ■ 

h3 



2(7 



n 



-{n + xf+ (l 



2(j 



n 



if j > i. Dividing this by (^) yields 

2(7 " * 



n 



n 



+ 1 



2(j - i) 



n 



n 



which proves part (b) of the theorem. 



□ 



A striking feature of the second part of Theorem 15.11 is that, for a fixed t, the 
probility that i < j is an inversion depends only on the difference j — i. This fact 
is not evident from the definition of the random process, but nevertheless it can be 
proved without invoking Theorem 15.11 



Proposition 5.2. Fix n > 2 and t > 0, and let ir E A n -\ be a product of t 
random transpositions. Then, the probability that 7Tj > ttj is an inversion (where 
1 < i < j < n) is only dependent on j — i. 

Proof. Let ir +k be the permutation obtained by adding k "modulo n" to all entries 
of 7r, i.e. 

7Tj + k if 7Tj + k < n, 

iTi + k — n if 7Tj + k > n. 
Clearly, for any 1 < i < j < n, we have 

Prob(7Tj > -Kj) = Prob{ir+ k > 7r+ fc ) + Prob^, > n — k) — Prob(7Tj > n-k). 
Now, fix k with l<i<j<j + k<n. It suffices to prove that 
(3) Prob(7r i+fc > TT j+k ) = Prob(7Tj > ttj). 



TT 



+ k 
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By symmetry of the random process, we have Prob(7Tj > n — k) = Prob(-7Tj > n — k) 



and hence Prob(7Tj > Wj) = Prob(-7r^" fc > 7rf k ). Equation ([3]) now follows from the 



observation that Prob(7Tj +/ t > itj+k) = Prob(7r 4 +fc > vr^ fc ). □ 

5.2. A definition and a lemma. The situation becomes a little more delicate for 
Weyl groups of type B and D, but in a way that is very similar for these types. To 
be economical we gather the bit of reasoning that is common for type B and D in a 
definition and a lemma. 

Definition 5.3. For any positive integer n, let Z n := {(i,j) : i,j G [— n, n] \ 

{o}, i*i + m. 

Lemma 5.4. Let n be a positive integer and x a real number. For every integer 
t > 0, define a function : I n — > M by the recurrence relation 



'ir=-S+ E v ej+ E »S 

\j'Mi\ 



together with the initial condition vfj = sgn(j — i). Then these functions are given 
by 



(t) j -i-Bgaj + Bgai t ( j - i - sgnj + sgn A , 

v) '■ = (2n - 2 + x) 1 + sgn(j - i) x\ 

hJ n — 1 \ n — 1 J 



Proof. Let T> n denote the real vector space of all functions v : I n — > R such that 
Vj t i = —Vij and V-j-i = Vij for any G X„. (We write instead of the more 
common v(i,j) since we may view T> n as the space of 2n x 2n matrices that are 
antisymmetric about the main diagonal and symmetric about the antidiagonal, if 
we disregard the entries on the diagonals.) 

Let Q be the linear operator on T> n defined by 

for any v G T> n and G I n , 
We have 

(Q 2 v)i,j = (Q^i'j + (Q v )hi' 
K'l^lil 

\i'\¥=\j\ \|i"l/b'l b"l#N'l / b'VNI \l*"l#b'l 

= (2n - 2) ^ ^2 Vi 'j" ~ ( v jj" JrV -3j") 

(i",j')eZn b"M|i| 
(now using the symmetry and antisymmetry properties) 

= (2n - 2) ^ ^ + uiji, 

\N"l^b'l b'"I^KI 
= (2n-2)(Q«) iii . 
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Thus, we have showed that Q = (2n — 2)Q which means that 



(Q + xlfv 



Q tfr Q 



(2n-2 + x) t —^ — + x t , 
v J 2n-2 \ 2n 

where I is the identity operator on T> n . The lemma now follows from the observation 
that 

(Qw (0) )ij = ^2 sgn ^ ~~ *') + E s S n (j' - = 2 U ~ i ~ s S n J + s S ni ) • 

\S'\¥\i\ 

□ 

5.3. The hyperoctahedral groups. Recall that B n is isomorphic to the subgroup 
of permutations n of the set [—n, n] \ {0} such that 7r_j = — 7Tj for all 1 < i < n (see 
e.g. [U Ch. 8]). In this representation, the set of reflections is 

{(», i)H> -i) : 1 < 1*1 < 3 < n} U {(i, -i) : 1 < i < n} 

in cycle notation. The length of ir equals its number of B-inversions, which are pairs 
with i, j £ [— n,n] \ {0} and j > \i\ such that 7Tj > ttj. 

Theorem 5.5. For n > 1 and £ > 0, Zei 7r^ 5e a product oft random reflections in 
B n . Then the following holds. 

(a) The expected length of it® is 

o / n 2 n(n+l) / 2\* n(n-2)( 4 2 X * 



2 3 \ n J 6 \ n n 2 

(b) G [— n, n] \ {0} and j > \i\, the probability that nf > 7rj is 

1 j-i-l + sgaif 1 2_y / j - i - 1 + sgni 1\ / £, 2 



2 2(n-l) V n / V 2(n-l) 2y \ n n 2 

For 1 < i < n, the probability that 7T_] > irf^ is 

1 _ 1 ( 2Y 

2 2V n, 

Proof. Part (a) of the theorem follows from part (b) by the simple observation that 

Ei>|i|(3-0 = 2n(n 2 -l)A 

It remains to prove part (b). For every i,j £ [— n,n] \ {0} and t > 0, let 

pg := Prob(vrf ) < Trf ) and ug := n 2 'pg- 

The following recurrence relation holds. 

(t+l) 2 (*) , / (t) (*h , / (*) (t)\ 

+ E (4j - u $) + E ( n S' - u S) if bi + ^d 

( 4 ) li'l^lil li'Mlil 



We also have the symmetry property _ i = ufj . 
,(*) ._ „.(*) „,(*) 



Now, define t> 4 j := u- • • — Ua\ for i, j € [— n, n] \ {0}. Clearly, v inherits the 

(t) (t) 
v ■ ■ = — f • - . 

3,t hJ 



symmetry property from u and has also the antisymmetry property vfj = —vf] . The 
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above recurrence for u is valid also for v and using the symmetry and antisymmetry 
properties it can be written 

(5) li'Mil ' b''l/KI 

(t+l) / N (t) 

Clearly, this implies that vf - = (re 2 — 2n) t for 1 < i < re, and if j > \i\ Lemma (f5T 
yields 

-.<*> - J 1 1 ■ SK "V - 2re)' + f 1 - J ' 1 ' S " U/ ') [» 2 - l» + 2)' 



M n — 1 V n — 1 

After dividing this by n 2t we obtain 
1 - 2p^ = pjj - p l hi 

j - i - 1 + sgni ( ^ _ 2\* + A _ j ~ j ~ 1 +sgni \ / _ 4 + _2_ x * 



n — 1 V n J \ re — 1 7 V re ' re 2 

if j > |i| and 

1 - 2 "'" = ( 1 -^ 

for 1 < i < n]. This proves part (b) of the theorem. □ 

5.4. The groups of type D. Recall that D n is isomorphic to the subgroup of 
permutations ir of the set [— re, n] \ {0} such that 7r_j = — 7Tj for all 1 < i < re and 
there is an even number of i £ {1,2,... , re} such that 7Tj < (again see e.g. [TJ 
Ch. 8]). In this representation, the set of reflections is 

~j) '■ 1 < HI < j < «} 

in cycle notation. The length of tt equals its number of D-inversions, which are 
pairs with i,j € [— n,n] \ {0} and j > \i\ such that 7Tj > ttj. 

Theorem 5.6. For re > 1 and t > 0, Zei ir® 5e a product oft random reflections in 
D n . Then the following holds. 

(a) T/ie expected length of tt® is 

D n(n-l) re(2n-l) / 2\ 4 n(n - 2) / 4^ ' 
^W = — 2 6 ^"nj 6 \_ re 

(b) Ifi,j G [—re, re] \ {0} and j > |i|j ^ e probability that irf^ > ir^ is 
1 j — i — 1 + sgn i ^ 2\' /j — i — 1 + sgn i 1 \ / 4 x * 



2 2(re-l) V »V V 2(re-l) 2/ \ n 

Proof. Part (a) of the theorem follows from part (b) by the simple observation that 

E^>|i|(3-0 = 2n(" 2 -l)A 

It remains to prove part (b). For every t > and G Zn> let 

pg := Prob(vrf } < Trf ) and ug := (n 2 - re)*pg • 
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The following recurrence relation holds: 

t45 hl) = (» a -«)«8+(«53-«S)+(«^-*-«S) 

(6) + E («&-«8)+ E («8-«8) 

K'MKUil b'^KI-bl 
We also have the symmetry property u_- _ i = ufj . 

Now, define vfj := ufj — for € Z. Clearly, inherits the symme- 

try property from it^ and has also the antisymmetry property u^*] = — v\ The 
above recurrence for it is valid also for u and using the symmetry and antisymmetry 
properties it can be written 

(7) 4^ = (n 2 - 5n + 4)«g + E ^ + E 

|i'l#b'l b'l/NI 

By Lemma 15.41 we obtain, for j > |i|, 

it) j — i — 1 + sgn i . 9 N . / 7 — i — 1 + sgn i \ , 9 . . 

v i i = ~ r^—fa - 3n + 2 * + 1 - J - ^-2— {n 2 - 5n + 4 *. 

M n - 1 V n - 1 / 

Dividing by (n 2 — n)* yields 

1 — ^Pj,i = Pi,j ~ Pj,i 

j -i- 1 + sgni ^ _ 2\* / _ j - j - 1 + sgnA / _ 4 X * 



n — 1 \ n ) \ n — 1 / \ n / 

if j > This proves part (b) of the theorem. □ 



6. Conclusion and ideas for future research 

We can sum up the current state of knowledge by the following table, listing the 

7 ( 



groups W for which we have exact expressions for Epf (t). 





R = simple reflections 


R = all reflections 


f = 


length 


A n , him), h(oo) 


A n , B n , him), G(r, l,n) 


if = 


absolute length 


him), h(oo) 


An, B n , D n , him) 



Perhaps the two most striking things with this table are that B n is missing from 
the upper left square and that not even A n appears in the lower left square! 

As Troili has pointed out [7J, the problem of computing E s "j{t) can be reduced to 
a heat flow process similar to that of the symmetric group, but it is not clear how 
to analyze it further. 

Computing Eg^it) is equivalent to computing the expected number of cycles 
(in the cycle representation of the permutation) after applying t random adjacent 
transpositions to the identity permutation. For now, we have no idea how to attack 
this natural problem. Note that it is not enough to keep track of the conjugacy class 
of the permutation as Eriksen and Hultman did in the proof of Theorem 13.31 

It is also worth noting that Egf^it) and EW e ,(i) is well-defined for any finitely 

generated Coxeter system. (Indeed, Troili made an effort to compute Egjit) in [7J, 
but unfortunately with an erraneous heat flow process.) 
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In addition to the length and the absolute length, there is at least one more 
"length" function that comes naturally with any Coxeter system (W, S), namely the 
descent number defined by 

d(w) = #{ S £S : £(ws) < £(w)}. 

(In fact, the ordinary length function is also a kind of descent number by the equality 
t(w) = #{t£T : t(wt) < £(w)}, see e.g. Hi Corollary 1.4.5].) What is Efjt) and 
Erp d (t) for a (in the second case, finite) Coxeter system (W, S)7 

7. ACKNOWLEDGEMENT 

This work was performed at KTH in Stockholm and was supported by a grant from 
the Swedish Research Council (621-2009-6090). 

References 

[1] A. Bjorner and F. Brenti. Combinatorics of Coxeter Groups. Springer, 2005. 

[2] M. Bousquet-Melou. The expected number of inversions after n adjacent transpositions. Discrete 
Math. Theor. Comput. Scl, 12:65-88, 2010. 

[3] N. Eriksen. Expected number of inversions after a sequence of random adjacent transpositions 
— an exact expression. Discrete Math., 298:155-168, 2005. 

[4] N. Eriksen and A. Hultman. Estimating the expected reversal distance after a fixed number of 
reversals. Adv. in Appl. Math., 32:439-453, 2004. 

[5] N. Eriksen and A. Hultman. Expected reflection distance in g(r, l,n) after a fixed number of 
reflections. Ann. Comb., 9:21-33, 2005. 

[6] H. Eriksson, K. Eriksson, and J. Sjostrand. Expected inversion number after k adjacent trans- 
positions. In D. Krob, A. A. Mikhalev, and A. V. Mikhalev, editors, Proceedings of Formal 
Power Series and Algebraic Combinatorics, pages 677-685. Springer, Berlin, 2000. 

[7] E. Troili. Forvantade avstand i Coxetergrupper (Expected distances in Coxeter groups). Mas- 
ter's thesis, KTH, 2002. In Swedish. 

Department of Mathematics, Royal Institute of Technology, SE-100 44 Stockholm, 
Sweden 

E-mail address: jonass@kth.se 



